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Abstract 



We develop the Tree-Loop Duality Relation for two- and three-loop integrals with multiple 
identical propagators (multiple poles). This is the extension of the Duality Relation for single 
poles and multi-loop integrals derived in previous publications. We prove a generalization 
of the formula for single poles to multiple poles and we develop a strategy for dealing with 
higher-order pole integrals by reducing them to single pole integrals using Integration By 
Parts. 
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1 Introduction 



The last 15 years have seen extraordinary progress in the analytical and numerical computation of cross 
sections in the Standard Model at one- and two-loops and even higher orders. Once considered the bot- 
tleneck of numerical applications, higher- than-tree-level corrections exhibit high complexity, partly due 
to the enormous number of Feynman diagrams needed, sometimes numbered in hundreds or even thou- 
sands, for important cross sections. Important efforts have been devoted to developing efficient methods 
able to boost forward the calculational capability both at the multi-leg and the multi-loop frontier. To- 
day, 2 — >• 4 processes at next-to-leading order (NLO), either from Unitarity based methods, JBElEl, 
or from a more traditional Feynman diagrammatic approach, [4], are affordable and are even becoming 
standardized. There has also been a lot of progress concerning next-to-next-to leading order (NNLO) 
calculations (MEU!. 

The importance of higher order corrections cannot be overstated. While leading-order (LO) predic- 
tions of multi-particle processes at hadron colliders in perturbative Quantum Chromodynamics (pQCD) 
provide, in general, a rather poor description of experimental data, NLO is the first order at which nor- 
malizations, and in some cases, the shapes, of cross sections can be considered reliable [9]. NNLO, 
besides improving the determination of normalizations and shapes, is also generally accepted to provide 
the first serious estimate of the theoretical uncertainty in pQCD. Despite the relatively smaller coupling, 
electroweak (EW) radiative NLO corrections might also be sizable at the LHC, ifTOlfTTf . Precision the- 
oretical predictions for background and signal multi-particle hard scattering processes, in the SM and 
beyond, are mandatory for the phenomenological interpretation of experimental data, and thus to achieve 
a successful exploitation of the LHC physics programme. 

In |fl2ll . a novel method was developed for the calculation of multileg one-loop amplitudes. Called 
the Duality Theorem, it applies directly the Cauchy Residue Theorem to one-loop integrals. The result 
can be represented by single cuts to Feynman diagrams, integrated over a modified phase space. (Note 
also lfl~3l where the author uses retarded boundary conditions to obtain a duality between loop and tree 
graphs). The Duality Theorem was extended in lfi~5l beyond the one-loop level, to two- and three- 
loops and it was shown how it can be extended to an arbitrary number of loops. The main feature and 
advantage of this approach is that at any number of loops, an amplitude can be written as a sum of 
tree-level objects, obtained after making all possible cuts to the lines of a Feynman diagram, one cut 
per loop and integrated over a measure that closely resembles the phase space of the corresponding real 
corrections. This modified phase-space, raises the intriguing possibility that virtual and real corrections 
can be brought together under a common integral and treated with Monte-Carlo techniques at the same 
time. In these papers the Duality Theorem was developed for diagrams and their integrals, that do not 
possess identical propagators. This possibility does not appear at one-loop for a convenient choice of 
gauge lfl~2l . but it is always present for higher order corrections. Identical propagators possess higher 
than single poles and the Duality Theorem developed so far, which is based on the Cauchy formula for 
single poles, must be extended to accomodate for this new feature. 

The purpose of this work is to extend the Duality Theorem to graphs and integrals with multiple 
poles at two- and three-loops and to present a procedure, a strategy for dealing with higher poles in an 
amplitude calculation that retains the features and advantages of the Duality theorem as detailed in |fl~5ll . 
The paper is organized as follows: In Section 2 we recall the basic definitions and results concerning the 
Duality Theorem at one- and two-loops and beyond, for integrals with single poles. In Section 3 we 
derive in detail the extension of the Duality Theorem to integrals with double and multiple propagators, 
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which exhibit multiple poles in the complex plane. In section 4 we present an alternative method to deal 
with multiple poles at two- and three-loops which retains the basic advantages of the Duality Theorem. 
It is based on Integration By Parts relations that allow to rewrite integrals with multiple poles in terms of 
integrals involving only single poles. 



2 Duality relation at one- and two-loops 

A general one-loop iV-leg, scalar integral, such as the one shown in Fig.[Q is written as: 

r N 

LU( Pl , P2 ,...,p N )= / Y[G F ( qi ), (1) 

J b i=i 

where 

g f(,Q*) = — TT~n W 

qf - mf + tO 

is the Feynman propagator. The four-momenta of the external legs are denoted Pi , i = {1,2, . . . N}. All 
are taken as outgoing and ordered clockwise. The loop momentum is l\, which flows anti-clockwise. 
The momenta of the internal lines q it are defined as 

qi = 2i+Pi,i, i e a x = {1,2, . . .JV} , P i,i= P i + ■ ■ -+Pi ■ (3) 
Throughout this paper we use the short-hand notation 

5(q i ) = 2ni9(q i , )5(q^-m 2 )=27ri5 + (qf-m 2 ), J . = -i J • (4) 

where the subscript + of 5+ refers to the on-shell mode with positive definite energy, > 0. Hence, 
the phase-space integral of a physical particle with momentum q if i.e., an on-shell particle with positive- 
definite energy, qf = m 2 , g i0 > 0, reads: 



d d q 



6{qi)---. (5) 



It was shown in [fT2l [T5l that using the Cauchy residue theorem the one-loop integral can be written in 
the form: 

f ~ n 

L {1) ( P i, P 2,-.-, P n) = ~z2 II G i>(Qi]Qj) > (6) 

i J ^ j=i 



Gd ft, Qj ^2 _ m 2 _ i0 v ( q . _ q ^ > ( 7 ) 



where 



is the so-called dual propagator, as defined in Ref. |I12| . with rj a future-like vector, 

Vo > 0, r, 2 = T)^ > , (8) 
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Figure 1 : Momentum configuration of the one-loop N -point scalar integral. 

i.e., a d-dimensional vector that can be either light-like {rf = 0) or time-like {rf > 0) with positive 
definite energy 770. The result in Eq. ©, contrary to the Feynman-Tree theorem (FTT) [fT71[T8l . contains 
only single-cut integrals. Multiple-cut integrals, like those that appear in the FTT, are absent thanks 
to modifying the original +i0 prescription of the uncut Feynman propagators by the new prescription 
— iO rj(qj — qi), which is named the 'dual' iO prescription or, briefly, the rj prescription. The dual iO pre- 
scription arises from the fact that the original Feynman propagator Gp{qj) is evaluated at the complex 
value of the loop momentum l\, which is determined by the location of the pole at q\ — mf + iO = 0. 
The iO dependence of the pole of Gp{qi) modifies the iO dependence in the Feynman propagator Gp(qj), 
leading to the total dependence as given by the dual iO prescription. The presence of the vector 77^ is 
a consequence of using the residue theorem and the fact that the residues at each of the poles are not 
Lorentz-invariant quantities. The Lorentz-invariance of the loop integral is recovered after summing 
over all the residues. Furthermore, in the one-loop case, the momentum difference r](qj — qi) is indepen- 
dent of the integration momentum £ x , and only depends on the momenta of the external legs (cf. Eq. ©). 

The extension of the Duality Theorem to two loops has been discussed in detail in [fT51[T6l . Here we 
recall the basic points. We extend the definition of propagators of single momenta to combinations of 
propagators of sets of internal momenta. Let ctk be any set of internal momenta qj with i, j e oi^. We 
then define Feynman and dual propagator functions of this set ctk in the following way: 

GVOfc) = Yl G p{qi) , G D {a k ) = ^2 5 (ft) ] I G D(qu qj) ■ (9) 

i&a k iea k jga fc 

By definition, Goictk) = S(qi), when ctk = {i} and thus consists of a single four momentum. At 
one-loop order, ctk is naturally given by all internal momenta of the diagram which depend on the single 
integration loop momentum t\, ctk = {1,2,..., iV}. However, let us stress that ctk can in principle be 
any set of internal momenta. At higher order loops, e.g., several integration loop momenta are needed, 
and we can define several loop lines ctk to label all the internal momenta (cf. Eq. (fT9l) ) where Eq. © 
will be used for these loop lines or unifications of these. We also define: 

G D (-a k ) = II °D(-qu -qj) , (10) 

where the sign in front of ctk indicates that we have reversed the momentum flow of all the internal lines 
in ctk- F° r Feynman propagators, moreover, Gp(— otk) = Gp(a k ). Using this notation the following 
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relation holds for any set of internal momenta 



G A {a k ) = G F (a k ) + G D (a k ) , (11) 
where G^(gj) is the advanced propagator: 

G A (qi) = ~ , (12) 

qf - mf - zO q ifi 

and 

GaM = 11 G A ( gi ) . (13) 

The proof of Eq. (fTTI) can be found in Ref. [15J. Note that individual terms in Gr>{(Xk) depend on the dual 
vector T], but the sum over all terms contributing to Gr>(oik) is independent of it. Another crucial relation 
for the following is given by a formula that allows to express the dual function of a set of momenta in 
terms of chosen subsets. Considering the following set 0n = «i U ... U a^, where 0n is the unification 
of various subsets ctj, we can obtain the relation: 

G D ( ai Ua 2 U ...Ua N ) II G D (a n ) JJ G F (a l2 ) . (14) 

The sum runs over all partitions of (3n into exactly two blocks (3jp and fi^ with elements « i; 
% G {1,...,N}, where, contrary to the usual definition, we include the case: 0^ = (3n, = 0- 
For the case of N = 2, e.g., where f3 2 = oli U a 2 , we have: 

G D (a 1 U a 2 ) = G D (ai) G D (a 2 ) + G D (ai) G F (a 2 ) + G F {a x ) G D (a 2 ) . (15) 

Naturally it holds that: 

N 

Gp^U^U ...Ua N ) = Y[G F (a l ) . (16) 

i=l 

Since in general relation (fT4)) holds for any set of basic elements a, L which are sets of internal momenta, 
one can look at these expressions in different ways, depending on the given sets and subsets considered. 
If we define, for example, the basic subsets a, to be given by single momenta g, , and since in that case 
Goifli) = 5 (qi), Eq. (fT4l) then denotes a sum over all possible differing m-tuple cuts for the momenta in 
the set /3n, while the uncut propagators are Feynman propagators. These cuts start from single cuts up 
to the maximal number of cuts given by the term where all the propagators of the considered set are cut. 
Using this notation, the Duality Theorem at one-loop can be written in the compact form: 

LV>(p 1 ,p 2 ,...,p N ) = - [ G D {a x ) , (17) 



where ai as in Eq. © labels all internal momenta qi . In this way, we directly obtain the duality relation 
between one-loop integrals and single-cut phase-space integrals and hence Eq. (fTTI) can also be inter- 
preted as the application of the Duality Theorem to the given set of momenta a x . It obviously agrees, at 
one loop, with Eq. ©. 

We now turn to the general two-loop master diagram, as presented in Fig. [21 Again, all external 
momenta are taken as outgoing, and we have p it j = Pi +Pi+i + ■ ■ - +Pj, with momentum conservation 
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Figure 2: Momentum configuration of the two-loop N -point scalar integral. 

Pi,n = 0. The label i of the external momenta is defined modulo N, i.e., pN+i = Pi- In the two-loop case, 
unlike at the one-loop order, the number of external momenta might differ from the number of internal 
momenta. The loop momenta are t\ and £ 2 , which flow anti-clockwise and clockwise respectively. The 
momenta of the internal lines are denoted by and are explicitly given by 

h+Pi,i ,ieati 
Qi= { h +Pi,i-i J e a 2 (18) 

k + h+Pi,i-i , i e a 3 , 

where oik, with k — 1, 2, 3, are defined as the set of lines, propagators respectively, related to the mo- 
menta qi, for the following ranges of i: 

ai = {0,l,...,r}, a 2 = {r + l,r + 2,. ..,/}, a 3 = + 1, 1 + 2, N} . (19) 

In the following, we will use for denoting a set of indices or the set of the corresponding internal 
momenta synonymously. Furthermore, we will refer to these lines often simply as the "loop lines". 



We shall now extend the duality theorem to the two-loop case, by applying Eq. (flTl) iteratively. We 
consider first, in the most general form, a set of several loop lines ct\ to a at depending on the same 
integration momentum i it and find 



/ Gp{ai U a 2 U ... U ctAr) — — Gd{&i U a 2 U ... U a^) , 

hi hi 



(20) 



which states the application of the duality theorem, Eq. (fTTT) . to the set of loop lines belonging to the 
same loop. Eq. (|20l) is the generalization of the Duality Theorem found at one-loop to a single loop 
of a multi-loop diagram. Each subsequent application of the Duality Theorem to another loop of the 
same diagram will introduce an extra single cut, and by applying the Duality Theorem as many times as 
the number of loops, a given multi-loop diagram will be opened to a tree-level diagram. The Duality 
Theorem, Eq. (|20l ), however, applies only to Feynman propagators, and a subset of the loop lines whose 
propagators are transformed into dual propagators by the application of the Duality Theorem to the first 
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loop might also be part of the next loop (cf., e.g., the "middle" line belonging to a 3 in Fig. [2]). The dual 
function of the unification of several subsets can be expressed in terms of dual and Feynman functions 
of the individual subsets by using Eq. (fl4l) (or Eq. (fT5T) ) though, and we will use these expressions to 
transform part of the dual propagators into Feynman propagators, in order to apply the Duality Theorem 
to the second loop. Therefore, applying Eq. (|20l) to the loop with loop momentum £i, reexpressing the 
result via Eq. (fT51) in terms of dual and Feynman propagators and applying Eq. (1201) to the second loop 
with momentum £ 2 , we obtain the duality relation at two loops in the form: 

L {2) (pi,P2,---,Pn) 

{-G D (ai) G F (a 2 ) G D (a 3 ) + G D (ai) G D (a 2 U a 3 ) + G D (a 3 ) Gd(-oii U a 2 )} . 

(21) 

This is the dual representation of the two-loop scalar integral as a function of double-cut integrals only, 
since all the terms of the integrand in Eq. ((2T]) contain exactly two dual functions as defined in Eq. ©. 
The integrand in Eq. ((21]) can then be reinterpreted as the sum over tree-level diagrams integrated over a 
two-body phase-space. 

The integrand in Eq. (|2TT) . however, contains several dual functions of two different loop lines, and 
hence dual propagators whose dual iO prescription might still depend on the integration momenta. This 
is the case for dual propagators Goiqi] qj) where each of the momenta qi and qj belong to different 
loop lines. If both momenta belong to the same loop line the dependence on the integration momenta 
in r](qj — q^ obviously cancels, and the complex dual prescription is determined by external momenta 
only. The dual prescription i](qj — q^) can thus, in some cases, change sign within the integration volume, 
therefore moving up or down the position of the poles in the complex plane. To avoid this, we should 
reexpress the dual representation of the two-loop scalar integral in Eq. (|2TI) in terms of dual functions of 
single loop lines. This transformation was unnecessary at one-loop because at the lowest order all the 
internal momenta depend on the same integration loop momenta; in other words, there is only a single 
loop line. 

Inserting Eq. (fT5T) in Eq. (12D) and reordering some terms, we arrive at the following representation of 
the two-loop scalar integral 

L (2) (pi,P2,- ■ -,Pn) 

= {G D {a 1 )G D {a 2 )G F (a 3 ) + G D (-ai)G F (a 2 )G D (a 3 ) + G*(a 1 )G D (a 2 )G D (a 3 )} , 

Jii Je.2 

(22) 

where 

G*(ai) = G F (a x ) + G D {a x ) + . (23) 

In Eq. (l22l) . the iO prescription of all the dual propagators depends on external momenta only. Through 
Eq. (1231) . however, Eq. ([221) contains also triple cuts, given by the contributions with three Goictk)- The 
triple cuts are such that they split the two-loop diagram into two disconnected tree-level diagrams. By 
definition, however, the triple cuts are such that there is no more than one cut per loop line a k . Since 
there is only one loop line at one-loop, it is also clear why we did not generate disconnected graphs at 
this loop order. For a higher number of loops, we expect to find at least the same number of cuts as the 
number of loops, and topology-dependent disconnected tree diagrams built by cutting up to all the loop 
lines afe. These results can be generalized at three-loops and beyond without any additional effort. The 
reader is refered to [fT5l for further details. 
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Figure 3: The two-loop N -point scalar integral with a double pole 

3 Duality Relation for multiple poles 

In the previous section we applied the Residue Theorem to one- and two-loop graphs that contain only 
single poles, i.e. no identical propagators. At one-loop this is always the case for a convenient choice 
of gauge [12J. However, at higher loops there exists the possibility of identical propagators, i.e. higher 
order poles. Obviously, we need to generalize the Duality Theorem to accomodate for such graphs. 
The first occurence of higher order poles is at the two-loop level, with the sole double pole generic 
graph shown in Fig. [3] The Duality Theorem can be derived for such graphs as well, using the Residue 
Theorem for multiple poles 



Res {2=2o} /(2;) 



(Jfe-1)! \dz k 



d 



k-l 



^) k f{z] 



(24) 



The derivation follows similar steps as with the single pole case and is independent of any particular 
coordinate system. We will derive an expression both in cartesian and light-cone coordinates, to demon- 
strate this independence. We start with the cartesian system. We write the Feynman propagator in a form 
that makes the poles explicit, i.e, 



G F (qi 



(?»o - 4 +) )(%o + Qio ] ] 



(25) 



where gjg = >/q^ + mf — iO is the position of the pole. Then, applying the Residue Theorem by 
selecting poles with negative imaginary part, we have 



R QS {Imq i0 <0}G F (qi 



N +) ) 3 



^ ldqo ^r 5+iqf - ml) - 



(26) 



The imaginary component of the new denominator 1/2 (q^ ) 2 is irrelevant, because it is always a positive 
quantity. We refer the reader to |[T2l where the calculation for the case of simple poles is explained in 
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more detail. Then, we assume the following Lorentz-invariant prescription of the residue 



ReS{ Img . 0<0 }Gp(gi 



dq - 



rj 



; 8 + {q 2 -m 2 ) 



(27) 



where r]^ = (r) , 0) is a future-like vector, i] > 0, in cartesian coordinates. Contrary to the one loop 
case, where numerators depending on the loop momentum do not modify the duality prescription, in the 
two loop and higher orders cases the derivative in the residue calculation introduced by the higher order 
poles act on every single term in the numerator and also on the remaining propagators. Let N(a k ) be a 
function of a set of momenta qi, with I £ a k . Then the residue of a double pole is given by 





Res { i mte)<0} \ G 2 F ( qi ) [ Y[G F ( qj ) ) N(a k ) \ = — - — ^ ( Y[G F ( qj ) ) N(a k ) 



Yl G D(qhQj] 



no=q ( i o ) 



N +) ) 2 



1 d 
[T) ~ 5Z( 2 ^o)G' D (g i ; qj) + — 



„(+) 



N{a k ) 



(28) 



which can be written as 



Res 



{lmg i0 <0} 



G 2 F { qi ) {j[G F { qj yj iVK)| = J dq 5 + (q 2 - m 2 ) (j[G D ( qi ; qj ) 



E— G D ( qi ; qj ) + -± 



d 



2r] qi dri qi 



N(a k ) . (29) 



In light-cone coordinates we choose our coordinates such that in the plus component complex plane 
the poles with negative imaginary part are located at: 

(+) _ Qi± + m l - ^ 



In these light cone coordinates the Feynman propagator reads 

Gf(Qi) 



with gj_ > . 



1 



and thus 



(30) 



(3D 



(32) 



Res{im qw <o}9(qi-)G F ( qi ) = , 

which, at first sight, seems to contradict equation Eq. (1271) . This contradiction can be resolved by taking 
into account the fact that in light cone coordinates, the dual vector t] is lightlike and therefore rj 2 = 0. 
Hence equation Eq. (1271 ) remains valid. Now, we are ready to calculate the residue of a double pole in 
light cone coordinates: 



Res 



{Im(jio<0} 



I 



dq+S+{q't ~ mf) | JJ G D (qi] qj) 



J d_ 

2r]q { drjqi 



N(a k ) 



(33) 
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where now r]^ = (rj + , r\- = 0, OjJ. Eq. (|33l ) has the same functional form as in Eq. (|29l) , although with 
a different dual vector 77. Thus we can generalize Eq. (1291 and Eq. (1331) to an arbitrary coordinate system 
and combining simple and double poles in a single formula we get in cartesian coordinates: 



G 2 F (q i )[l[G F (q j )\N(a k ) = 



2 (TO 



j d_ 



(34) 



Equation (|34l) . is the main result of this section. It extends the Duality Theorem to integrals with 
identical propagators or, to put it differently, with double poles in the complex plane. For the case of 
the generic two-loop graph in Fig. [31 this result can be seen as an extension of Eq. (fT5l) . If we have 
two groups of momenta, a k , a 2 , one of which contains the double propagator, i.e. a k = {q n = £1 + £2} 
and a 2 = {q 2 = £2, qz = h + Pi, ■ ■ ■ , Qn-i = h+ Pi,n-i, <?2 = 4}, and we denote by a' 2 a group that 
contains all the momenta of a 2 leading to single poles, namely 



«0 



{q 2 = £ 2 , q 3 = £ 2 +Pi, . . ■ , q n -i =£2+ Pi,n-i}, then we can write: 



G D {a k Ua 2 ) = 5{q 2 ) G D {q 2 ;qj) 



2(7752 



9 n 

.rr 1 vq2 



+ ~ 5 ^)G 2 D (q t ;q2) ( Yl G D ( qi ;q 3 ) 



(35) 



This result states that for the case of a double pole, one follows the usual procedure of cutting every 
propagator line once, including the double propagator, and transforming the rest of the propagators to 
dual propagators. A similar formula can be derived for the case of multiple (triple and higher) poles. 
The calculation of the residue of a multiple pole introduces, however, contributions with powers of dual 
propagators. In absence of a general transformation formula analogous to Eq. (fT4l) . it is not possible 
to rewrite Eq. (|34l) in terms of dual propagators whose dual +2O prescription depends on the external 
momenta only. For that reason, we will present in the next section a different strategy for dealing with 
higher order poles based on the reduction of the integral using Integration By Parts. 



4 Reducing to single poles with IBPs 

In this section, we discuss a different approach to the generalization of the Duality Theorem to higher 
order poles. We will use Integration By Parts (IBP) [120112711 to reduce the integrals with multiple poles to 
ones with simple poles. We emphasize the fact the we do not need to reduce the integrals to a particular 
integral basis. We just need to reduce them "enough", so that the higher order poles disappear. 

To give a short introduction to the method and establish our notation, let us consider a general m- 
loop scalar integral in d dimensions, with n denominators D 1 , . . . , D n raised to exponents a 1; . . . , a n and 
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external momenta pi, . . . 



,Pn- 




(36) 



If we notice that 




St, 



d f 



(37) 



ds» Dl 1 ■ ■ ■ D% 



0. 



where s M = , ■ ■ ■ , the integrand being a total derivative with respect to the loop momenta, we 
can find relations between scalar integrals with different exponents aj. This will allow us to express 
integrals with exponents larger than one, in terms of simpler ones. In effect, we will be able to write 
integrals with multiple poles in terms of sums of integrals with simple poles. In the numerator of the 
integrand of Eq. (1371) we can use t M = . . . , p^,. . . , p 1 ^, to obtain a system of equations that relate 
the various integrals. For simplicity, when refering to an IBP we will use the shorthand notation: 



to denote Eq. (1371) . The differentiation will raise or leave an exponent unchanged, while, contractions 
with the loop and external momenta in the numerator of the integrand, can be expressed in terms of 
the propagators to lower an exponent. Often times though, this is not possible, leaving scalar products 
of momenta, which cannot be expressed in terms of denominators. These are called Irreducible Scalar 
Products (ISP). We will consider ISPs as additional denominators, = ^ • pj , with a negative index 
dj. We use the notation: 



to denote a generic two-loop integral with n propagators raised to an arbitrary integer power, with 
Di = G~ F 1 {qi) = qf — + iO and q- L denotes any combination of external and loop momenta. In the 
following the prescription +i0 for the propagators is understood. We will use the symbol to denote 
the raising of the index a* by one i.e. 1 + F(ai, a 2 , • • • a n ) = F{a% + 1, a 2 , . . . , a n ) and the symbol 
to denote the lowering of the index a { by one i.e. 2~F(ai, a 2 , • • • a n ) = F(a\, a 2 — 1, . . . , a n ). A 
combination of the two means that the operators apply at the same time i.e. l + 2~F(a!, a 2 , • • • a n ) = 
F(a\ + 1, a 2 — 1, . . . , a n ). In the following we will use two automated codes, for the reduction, FIRE 
Q21 , a MATHEMAT I CA package for the reduction of integrals andREDUZE 2 Q , a package written 
in C++, using GiNaC 

4.1 The case for two-loop diagrams 

The only generic two-loop scalar graph with iV-legs and a double propagator is shown in Fig. [3l The 
simplest case is the two-point function with massless internal lines. The denominators are: 



* Since the most obvious first approach seems to be to try to express the integrals with multiple poles in terms of the same 
integrals with only single poles, c.f. Eq. (HTI) . we used, in addition to the "usual" version of REDUZE 2, in some cases a 
special patch for REDUZE 2 which provides a modification of its integral ordering in the final result. This modified version 
of REDUZE 2 delivered the results for the integrals in this desired form stated in the subsequent sections, while we used the 
normal version of the integral ordering for the remaining cases. Note that we also calculated explicitely the relations obtained 
from the modified version, in the easiest cases of the massless two- and three-loop integrals which can be built by insertion 
of the the massless one-loop two-point function, and found agreement. 




(39) 



= t\ , D 2 = ej , D 3 = (£ 2 +p) 2 , J D 4 = (£ 1 + £ 2 ) 2 , D 5 =h-p, 
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where we have introduced an ISP as an additional denominator Q. In our notation, the integral we want 
to reduce is F(12110) and to this end we use the six total derivatives 



dL 1 ' dt; 



, • P , i,J = 1,2. (40) 



Applying these IBPs on ^(0102030405) we get a system of recursive equations. Using specific values 
for the exponents Oj we can solve this system and obtain F(12110). For this particular case, we solve 
the system explicitely and the reader is referred to the appendix [A] for details. Finally we arrive at: 

F(12110) = F(llllO) , (41) 

(l + e)a 

where s = p 2 + iO, a result which contains only single poles and can be treated using the Duality 
Theorem [15]. For the rest of the cases below and in the three-loop case in the next section, we have 
used FIRE and REDUZE 2 to perform the reductions and check our results. For three external legs 
Pi = p\ — 0;Pl — (Pi + P2) 2 and massless internal lines, we have the denominators: 



D 1 = £i , D 2 = q , D 3 =(£ 2 + Pl ) 2 , D 4 =(£ 2 + Pl +p 2 y , D 5 = {h+£ 2 y , 

A> = £i-pi , D 7 = £ x ■ p 2 , 

where, the last two are the ISPs that appear in this case. The integral we want to reduce is F(1211100). 
We use eight IBPs: 

W • iJ = 1 ' 2 - (42) 

A similar analysis to the one above, gives: 

F(1211100) = , 36 r F(lllllOO) = - 3(1 ~ 3e)(2 ~ 3e) F(1001100) , (43) 
(l + e)s e(l + e)s 3 v ; 

where s = (pi + p 2 ) 2 + iO, which, again contains only single poles and can be treated with the Duality 
Theorem. 

The inclusion of masses does not affect the general picture of the reduction. It solely introduces 
numerators in some integrals after the reduction is done. But, as we have stressed already, the application 
of the Duality Theorem is not affected by numerators since it only operates on denominators [fT5l . As 
an illustrative example, let us consider the two-loop graph with two external legs and one massive loop 
(see Fig. [3]). For the case of the left loop being massive (related to l 2 ), with mass m, the denominators 
involved are 

D 1 = l\ , D 2 = l\-m 2 , D 3 = (£ 2 + p) 2 - m 2 , D A = (£, + l 2 ) 2 - m 2 , 
with the addition of the irreducible scalar product 

D 5 = £t-p, (44) 



^For REDUZE 2 the corresponding propagator is added and used as input instead. 
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needed to perform the reduction. Using the same IBPs of Eq. (1401 ), the result of the reduction, with FIRE 
is: 

F(12110) = (( - l) ^ + 2m ' (9e - - 3)S + im \ 3 + 2f)( ~ 1 + 2f)l F(00110) 

2e(2e - l)m 4 s (4m 2 - s) 2 

-J- 2 F ( 10110 ) - TTZ^7 F ( 1 - 1U0 ) 



2e(2e-l)m 4 s(4m 2 -s) 2 2(2e-l)mV 

(e - l) 2 (2m 2 - s) . . (e - 1) (4em 2 + 2m 2 -s) , 

-F(01010) - \. n M 2 r^(OlllO) 



(2e - l)m 4 s (4m 2 - s) v 7 2(2e - l)m 4 (4m 2 - s) 
2(e _ 1)K _ s)(W _ 6 _ 3 ^) 

e(2e - l)m 4 s (4m 2 - s) 2 

with s = + iO and 

A = (1 - e) [es + 2(3 - 8e)m 2 ] s 2 + 2(1 - 2e)m 4 [2(3 - 4e)m 2 - (6 - 5e)s] . (46) 
The reduction generates two integrals with a numerator, namely 



F(l - 1110) 
F(1011 - 1) 




1)2 ^2 



■2 



m 




D1D3D4 ' 
ii-p 



_ D,D 3 D 4 

but the double poles have now disappeared. The result with REDUZE 2 reads: 

(e - 1) (4em 2 + 2m 2 - s) . 

F 12110 = — - — - — — — r — —-j- — -^F OHIO 

v ' 2 2e - 1 m 4 4m 2 - s) V ' 



3(e - 1) (8em 4 - 12em 2 s + es 2 - 4m 4 + 4m 2 s) 

+ ^ r (1 — 1110) 

2e(2e - l)m 4 s (4m 2 - sf 



+ ^F(lOllO) 

2e(2e - l)m A s (Am 2 - s) 2 

(e — 1) (8e 2 m 2 s — 2e 2 s 2 — 16em 4 + 6em 2 s + es 2 + 12m 4 — 6m 2 s) , 

+ 5 -F(01010) 

2e(2e - l)m 4 s (4m 2 - s) 



(47) 



where A is given by Eq. (|46l ). Despite the appearence of different integrals the two results are of course 
equivalent. This is because, the integrals F(00110) and F(01010), in the result obtained with FIRE, are 
identical (as can be seen by shifting the loop momenta), so the sum of their coefficients gives exactly 
the coefficient of the result obtained with REDUZE 2. The same argument applies for the integrals 
F(1011 — 1) and F(l — 1110). The appearence of the numerators does not affect the application of the 
Duality Theorem for integrals with single poles as was detailed in [fT5l . For the case of the right loop in 
Fig.[3]being massive (related to £1), we have the denominators: 

D x =i\-m 2 , D 2 = i\ , D 3 = (£ 2 +p) 2 , D 4 = (£ 1 +£ 2 ) 2 - m 2 , D 5 = i x ■ p . 
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Using the IBPs from Eq. (|40l , we get with FIRE: 

(32e 2 m 4 + 8e 2 m 2 s + e 2 s 2 - 32em 4 - llem 2 s - es 2 + 6m 4 + 3m 2 s) 
F(12110) = ^ F(10110) 

(e - 1) (16e 3 m 2 + 4e 3 s - 20e 2 m 2 - 8em 2 - 7es + 3m 2 + 3s) 

1 (10010) 



6(2e- l)(2e+ l)m 4 s 2 
e - 1)f F(1011-l)-(i^l^F(01110) 



3m A s ' 2m 2 s 
(e - 1) (12em 2 + es - 3m 2 



6m 4 s 2 

(e — 1) (6em 2 + es — 3m 2 ) 



6m 4 s 

and with REDUZE 2: 



^F(l - 1110) 
F(11100) , (48) 



F(12U0) = _M(^ f|1 . mo) 

64e 2 m 4 + 16e 2 m 2 s + e 2 s 2 — QAem 4 — 22em 2 s — es 2 + 12m 4 + 6m 2 s , 

+ 12m^ F(1 ° 110) 

(^-l)(12em 2 + es -6r^) F(llloo) 



6m 4 s 

(e - l)(2e - 3) (16e 2 m 2 + 2e 2 s + 4em 2 + 3es - 2m 2 - 2s) 
12(2e- l)(2e+ l)m 4 s 2 



F(IOOIO) . (49) 



For the case of the double pole, two-loop graph, with three external legs and one massive loop, we have 
the denominators: 

D 1 = i\ , D 2 = t\-m 2 , D 3 = (£ 2+Pl ) 2 -m 2 , D 4 = (£ 2+Pl +p 2 ) 2 - m 2 , D 5 = (£ 1+ £ 2 f-m 2 , 

D 6 = £ 1 -p 1 , D 7 = £ x ■ p 2 . 
Using the IBPs from Eq. (021) we get with FIRE: 

(e-l){(l + 4e)s 2 -4em 2 s(ll-2e)-8m 4 (4e 2 -8e-l)} . 

F 1211100 = ^ Jxy — '- i '- =-i — F(OOOllOO) 

8e(2e - l)m 6 s (Am 2 - s) 2 

A 

— ^F(lOOllOO) 

2e(2e - l)m 6 s (Am 2 - s) 

(8e 3 -12e 2 + 4e-l)(e-l) . 

+ , e M -F 0010100 

8e(2e - l) 2 m 6 s v ; 

{( ~ ' ' [F(OllllOO) - F(lOlllOO)] + ^rf-^r c -F(OlOllOO) 



2(2e - l)m 4 m 2 s (Am 2 - s) 

( £ _ 1) ( 2em 2 es m 2) 

+ 26( 2e - l)m*s 2 ^i- 101100 ) 

+ ie - in8em '- 2 r- 2 6m \ +S) F(0100l00) 
2(2e - l)m 6 s (Am 2 - s) v ; 

- * - V ^ ~ «> ^ + m ' a & - 1) + ^ ~ ^FdOOll - 10) , (50) 
e(2e - l)m 6 s 2 (Am 2 - s) 2 
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where s = (pi + P2) 2 + *0, and: 

A l = e (l - e)s 3 + m 2 (-l + e)(-l + 9e)s 2 + m 4 (l - e)(5 + 6e)s + 2m 6 (l + 2e)(-3 + 4e) , (51) 
while, with REDUZE 2, we get: 

F(12moo) = - 2(2^TR F ( 01111tx) ) 
+ 2(27fTR F(1011100) 

m z s (4m z — s) 

3(e - 1) (4em 4 - 8em 2 s + es 2 + 2m 4 + m 2 s) 

— ^ b (1 — 1U11UU) 

2e(2e - l)m 6 s (4m 2 - sf 

- — ? F(1001100) 

2e(2e - l)m 6 s (4m 2 - sf 

A* 

+ ^F(OIOOIOO) , (52) 

4e(2e - l) 2 m 6 s (4m 2 - s) 2 v ; ' 



where 



and 



A 2 = 16e 2 m 6 - 6e 2 m 4 s + 9e 2 mV - e 2 s 3 

— 4em 6 + em 4 s — 10em 2 s 2 + es 3 — 6m 6 + 5m 4 s + m 2 s 2 , (53) 



A 3 = (e - 1) (I28e 4 m 4 - 64e 4 m 2 s + 8e 4 s 2 - 256e 3 m 4 + 112e 3 m 2 s - 12e 3 s 2 + 192e 2 m 4 

- 92e 2 m 2 s + 8e 2 s 2 - 40em 4 + 26em 2 s - es 2 - 12m 4 + 4m 2 s - s 2 ) . (54) 



The cases with additional external legs can be treated in a similar manner. It can always be reduced to 
sums of integrals with single propagators at the expence of introducing numerators. Although no formal 
proof exists, in all cases studied so far it has been possible to reduce to integrals where only single poles 
appear. The generality of this result seems plausible [fT4l . 

Our strategy is now clear. For a two-loop calculation, first we reduce all double pole graphs using 
IBPs or any other method. The remaining integrals all contain single poles and can be treated using the 
Duality Theorem at two-loops. The appearence of vector or tensor integrals does not spoil this strategy 
since the Duality Theorem for single poles, affects only the denominators of the integrands. 



4.2 The case for three-loop diagrams 

For three-loop graphs there exists one topology with a triple propagator and a number of topologies with 
a double propagator. All topologies are shown in Fig. |H The arguments for the two-loop case are valid 
here as well. We first reduce the multiple pole integrands by using IBPs until we have integrals with only 
single poles (possibly with numerators) and then we can then apply single-pole Duality Theorem as it 
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M) to (/) 




Figure 4: Master topologies of three-loop scalar integrals with multiple powers of internal propagators. 
Each internal line already dressed with some external leg can be dressed with additional external lines. 



was described for the three-loop case in Ref. [15J. In the following, we show explicitly the reduction of 
the two-point function for the different topologies and for massless internal lines. We use the notation: 

F( ai a 2 ---a n )= } (55) 

to denote a generic three-loop integral with n propagators raised to an arbitrary integer power, with 
Di = Gp 1 (q i ) = q 2 — mf + iO or Di equal to any ISP and qi any combination of external and loop 
momenta. For the rest of this section the prescription +i0 for the propagators is understood. We also 
have s = p 2 + iO. The IBPs to be used for the reduction are: 

i,i = 1,2,3 (56) 

In the following, we present first the result obtained with REDUZE 2 and then with F IRE. For the single 
triple pole graph (a) in Fig. ©, we have the following expressions: 

(a) The denominators used are: 

D 1 = t\ , D 2 = i\ , D 3 = l\ , D 4 = (£ 2 - P ) 2 , D 5 = (e x - £ 2 ) 2 , D 6 = {£ 3 - i 2 ) 2 , 
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D 7 = £ l -p, D 8 = £ 3 -p, D 9 = £ l -£ 3 , 

with the result: 

F(131111000) = 2e(-l + 4e) F ( mm0 oo) = ghl+Mzl+i^ F(lOllllOOO) . 
1 ; (l + e)(l + 2e)s 2 V ; (l + e)(l + 2e)s 3 V ; 

(57) 

For the graphs with doubles poles, (b)-(g), Fig. ©, we find: 

(b) The denominators are: 

D 1 = P l , D 2 =f 2 , D 3 = 4, D 4 = (£ 2 -p) 2 , D 5 = (£, - £ 2 f , D 6 = (£ 3 - £ 2 + pf , 

D 7 = i x -p, D 8 = £ 1 -£ 3 , D 9 = £ 3 -p, 

with the result: 

3(-l + 4e)(l + 3e) , 
F(121211000) = — - j^- 2 ^F(llllllOOO) 

= 6(-2 + 3 t )(-l + 3 f )(l + 3 f )(-3 + 4 f )(-l + 4 t ) Q101 

e 2 (l + e) 2 (-l + 2e)s 4 V ; 

(58) 

(c) The denominators are: 

A = £\ , 7J 2 = £\ , 7J 3 = ^ , D 4 = (£ 3 +p) 2 , D 5 = (£ 3 - £ 2 f , D 6 = {£, - £ 2 ) 2 , 

D 7 = l x ■ p , D 8 = £ 2 ■ p , D 9 = l x ■ £ 3 , 

with the result: 

2e(-l + 4e)(-l + 3e) , 
F(122111000) = ( \ + 2e)( ; l +e)2g2 V (llllllOOO) 

= 2( - 2 + 3e)( - 1 + g ( - 3 + 4 4 £)( " 1 + 4£) F(lOOlllOQQ) . (59) 
e(l + e) 2 (l + 2e)s 4 V ; 

(d) The denominators are: 

D 1 = £l, D 2 = £ 2 2 , D 3 = f 3 , D 4 = (£ 3 -p) 2 , D 5 = (£ 2 + £ 3 - £ x f , 
D^^-p, D 7 = £ 2 -p, D 8 = £,-£ 2 , D 9 = £ x -£ 3 , 

with the result: 

F(112110000) = (~ 1 + 2e ) F(lllllOOOO) = (~ 3 + 4e ) F(llOllOOOO) . (60) 

(e) The denominators are: 

D 1 = £j, D 2 = £l, D 3 = £j, D 4 = (£ 3 -p) 2 , D 5 = (£, + £ 2 ) 2 , D 6 = (£ 2 + £ 3 ) 2 , 

D 7 = i x -p, D 8 = £ 1 -£ 3 , D 9 = £ 2 -p, 

with the result: 

. (-l + 4e) , , (-2 + 3e)(-3 + 4e)(-l +4e) , 
F(112111000) = j- — pF(llllllOOO) = - '-j- — - 3 - F(lOOlllOOO) . 

(61) 
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(f) The denominators are: 

D 1 = ( 2 , D 2 = £ 2 2 , D 3 = f 3 , D, = (£ 2 +p) 2 , D 5 = {£, + £ 2 ) 2 , D e = + £ 3 ) 2 , 

D 7 = (e 3 -e 2 ) 2 , D % = i x . P , D 9 = e 3 . P , 

with the result: 

F(121111100) = ^ - F(lllllllOO) 

= 2(-2 + 3 e )(-l + 3A<dt0tl ± ^ [f (QonilOOO) + F(IQQIOIIOO)] 
2(-l + 2e) 2 (-l + 4e) , 

(g) The denominators are: 

A = ^ , D 2 = i\ , D 3 = t\ , D A = (£ 3 -p) 2 , 

D 5 = - £ 2 ) 2 , D 6 = (£ 3 - £ 2 ) 2 , D 7 = (£ 3 -l 2 - pf 

D 8 = (h~£ 3 ) 2 , D 9 = (h-p) 2 (63) 



with the result: 



F(i2imioo) = T^^0 mnni - 10) 



+ (iW^ F(iiiiiiioo) ■ (64) 

The difference between the two results is due to the fact that the second is expressed in terms of basis 
integrals while the first is expressed in terms of integrals with single poles of the same type as the 
multiple pole integral (in effect the first result can be further reduced to the second). Since we do not 
seek a particular basis for our reduction, as was stressed earlier, both results are equally useful as far as 
application of the Duality Theorem is concerned. 



5 Conclusions 



We have extended the Duality Theorem to two- and three-loop integrals with multiple poles. A Lorentz 
invariant expression for the residues of double poles has been derived, which can be extended straightfor- 
wardly to triple and, in general, multiple poles. In the absence of a systematic procedure to express the 
dual -MO prescription in terms of external momenta exclusively, as in the case of single poles, we have 
explored an alternative approach. We use IBP identities to reduce the integrals with identical propaga- 
tors to ones with only single poles. Therefore, the essential features of the loop-tree duality now remain 
intact. We reiterate that our goal is not to reduce everything to some set of master integrals. Rather, we 
reduce the integrals until there are no multiple poles left. Then, we can use the Duality Theorem in its 
original form for single pole propagators, to rewrite them as integrals of a tree-level object over a modi- 
fied phase-space. The appearence of additional tensor integrals, due to the reduction, does not affect our 
procedure, since applying the Duality Theorem in its single-pole version, only cuts propagators, leaving 
the numerators of the integrals unaffected. 
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A Proof of the reduction of Eq. (I4H) 



Here we solve the system of equations, explicitely, to arrive at Eq. (14TT) . We note that we are not aiming 
for a full reduction to a set of master integrals but rather to reduce the multiple poles to single poles. 
Therefore, any integral which has single propagators is to be considered known. 

Using the IBPs, Eq. (l40l) . on the generic integral F(aia 2 a 3 a 4 a 5 ), we get the system of equations: 

d — 2a% — a 4 — a 5 — a 4 4 + l~ + a 4 4 + 2~ = , (65) 

a 1 -a 4 + ^sa 5 5 + + a 4 4 + (l~ - 2~) + ail + (2~ - 4") + ^a 5 5 + (2~ - 3~) = , (66) 

a 4 4 + (s + 2 - 3 - 2 5") - sa 5 5+ - 2ail+5~ = , (67) 

a 2 - a 4 + a 2 2 + {l~ - 4") + a 3 3 + (l" + 2 - 4 - 2 5") + a 4 4+(2~ - IT) = , (68) 

d - 2a 2 - a 3 - a 4 + a 3 3 + (s - 2") + a 4 4 + (l~ - 2") = , (69) 

a 2 - a 3 + s(a 2 2+ - a 3 3 + + a 4 4 + ) + a 4 4 + (2~ - 3" - 2 5") + a 3 3+2~ - a 2 2+3~ = , (70) 

where s = p 2 + i0. The appearence of the operator 5 signals that we have the ISP i\ -p in the numerator 
of an integral. As long as these integrals possess single propagators, we will not reduce them further but 
consider them known. We also note that a lot of the integrals that appear after setting particular values 
to the parameters a, in this system, are zero in dimensional regularization (in the massless case). Let us 
start by setting a 2 = 2, a± = a 3 = a 4 = 1, a 5 = —1 in (|66l ). We get: 

F(2111 - 1) - F(1112 - 1) - -sF(12110) - ^F(llllO) = . (71) 



Taking the sum of (1691 ) and (1701) and setting a\ — 1, a 2 — a 3 — a 4 = 1, a 5 = we get: 

(d - 4)F(11110) + sF(12110) + sF(11120) - 2F(1112 - 1) = . (72) 
Taking the difference between (l69l) and ([70]) and setting a\ — 2, a 2 = a 3 = a 4 = 1, o 5 = we get: 

- 2F(2111 - 1) - sF(12110) + sF(11210) - F(10210) = . (73) 
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Finally, setting a 1 = 2, a 2 = a 3 = a 4 = 1, a 5 = —1 in (1651 ), we get: 

(d - 4)F(2111 - 1) - F(1112 - 1) + F(2012 - 1) = . (74) 

The integrals F(10120) and F(11210) and F(2012 — 1), in this system of equations, can be computed 
simply by taking further, appropriate combinations of Eqs. (l63T)- (l70l . Setting cti = a 2 = a 3 = a 4 = 
1, a 5 = in d65j we get F(10120) = (3 - d)F(llllO). It also holds that F(10120) = F(10210), 
as can be seen by making the shifts in the loop momenta i 2 — > — ^2 — V- From Eq. (|69l) . by 
setting a\ = a 2 = a 3 = a 4 = 1, a 5 = and using the value of F(10120), we get: sF(11210) = 
(10 — 3g?)F(11110). Finally, adding Eqs. (1651) and (|66l) and setting a\ = a 2 = 1, 03 = 2, a 4 = 1, as = —1 
we get: F(2012 - 1) = (d- 3)(d-4)F(11110). The rest of the system of equations (f7Tb-d74l) can now 
be solved sequencially, arriving at: 

F(12110) = ( . 3 , d ~\ 0) F(11110) , (75) 
[d — b)s 

a result which, after putting d = 4 — 2e, agrees with Eq. (|4~TT) . 
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